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A gauge-invariant Wigner quasi-distribution function for charged particles in classical electromag- 
netic fields is derived in a rigorous way. Its relation to the axial gauge is discussed, as well as the 
relation between the kinetic and canonical momenta in the Wigner representation. Gauge-invariant 
quantum analogs of Hamilton- Jacobi and Boltzmann kinetic equations are formulated for arbitrary 
classical electromagnetic fields in terms of the 'slashed' derivatives and momenta, introduced for 
' this purpose. The kinetic meaning of these slashed quantities is discussed. We introduce gauge- 

invariant conditional moments and use them to derive a kinetic momentum continuity equation. 
£SJ ■ This equation provides us with a hydrodynamic representation for quantum transport processes and 

a definition of the 'collision force'. The hydrodynamic equation is applied for the rotation part 
of the electron motion. The theory is illustrated by its application in three examples. These are: 
Wigner quasi-distribution function and equations for an electron in a magnetic field and harmonic 
potential; Wigner quasi-distribution function for a charged particle in periodic systems using the kq 
qq ■ representation; two Wigner quasi-distribution functions for heavy-mass polaron in an electric field. 
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I. INTRODUCTION 

In 1932 Wigner introduced a possible method for describing a quantum system in terms of a phase space quasi- 
distribution function which is now known as the Wigner quasi-distribution function. It allows quantum mechanical 
expectation values to be written in the form of phase space integrals similarly to the corresponding expectation values 
in classical mechanics (see, e.g., Balazs and Jennings, 1984; Hillery et al, 1984). The Wigner quasi-distribution func- 
tion is closely connected with the Weyl (1931) correspondence principle (see also Schwinger, 1960), which has been 
Ch , originally formulated for systems free of external electromagnetic fields. 
The Wigner quasi-distribution function, 

^ \ Pw(Pcan, X) — pwi^can, Pcan, T, R) , 

* *~H i 

is defined as an expectation value of the so-called Wigner operator, pw{Pcan,X) obtained by a coordinate transfor- 
mation of a single particle density operator, p(x, x') = p(t, r, t' , r'). This transformation Wigner (1932) was chosen in 
d ' the form 

_ _ i 

Pw (Pcan,X)= [ &*ye-W p ^p{y,X) (1) 



where X — (x + x')/2, y — (x' — x) and the subscript can indicates that P in this expression is the variable conjugated 
to the coordinate X. Although the transformation ([l]) seems quite obvious, it has a highly significant drawback of 
producing not gauge-invariant quantities for charged particles in an electromagnetic field. The resulting conjugate 
variable P can — {^cam Pcan) does not necessarily have a meaning of the kinetic four momentum. 

Here lies an important difference between the classical and quantum descriptions. In classical mechanics we deal with 
the well-defined coordinate and momentum of the particle. The electromagnetic field is introduced in the Hamilton 
function of the system simply by considering the canonical momentum 

e 

Pcan Pkin ~\~ ~ ^(^) 
C 
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instead of the kinetic momentum. Here A^(x) = (<p(x), A(x)) is the four vector-potential of the applied electro- 
magnetic field. While it is well known that no problems with defining gauge-invariant quantities arise in classical 
mechanics, the uncertainty principle does not allow one to define simultaneously a coordinate and a momentum of a 
quantum particle. Quantum mechanics must, therefore, deal with a two-point density operator p(x, x') — p(t, r, t', r'). 
Then a gauge transformation converts the potential A(x) differently in these two different points, making an intro- 
duction of the canonical momentum not so harmless and even ambiguous procedure. The result is that equation ([!]) is 
not the only, nor is it the best definition of the Wigner quasi-distribution function for a system in an electromagnetic 
field. 

Nonetheless, the transformation (nl) has been used many times in the past for systems in electromagnetic fields. 
This may lead to certain difficulties^] in calculations and in interpretation of the results (see, e.g., Mahan, 1987; 
Mahan, 1990; Rammer, 1986, and the discussion in Levanda and Fleurov, 1994) since terms explicitly depending on 
the potential may appear, e.g., in kinetic equations. 

Fleurov and Kozlov (1978) (see also Al'tshuler, 1978) applied a transformation of the form 

p w (P,X) = [d*ye-* v "( p »-° A ^ x »p(y,X), (2) 



which produces explicitly gauge-invariant quantities and equations in the linear approximation or in a nonlinear case 
of a static homogeneous electric field (Tugushev and Fleurov, 1983; Serimaa et al, 1986 ). 

The first who tried formulating a gauge-invariant Wigner quasi-distribution function for an electron in classical 
electromagnetic fields was Irving (1965) who proposed to employ the transformation 

pw(P,X) = 

J dVxp f-ipy + ii J y^A^X + sy)d)j p(y,X) (3) 

which is truly gauge- invariant. It is worthwhile to pay a particular attention to the phase factor in ([|) emerging from 
the proper time method of Schwinger (1951). This phase factor is related to the dynamics of the charged particle in 
the electromagnetic field. Really, the phase factor of an electron wavefunction which transforms as 



under the local gauge transformation 



i)(x) -> iP(x)e iie / hc ^W (4) 



A„(x) -> An(x) - M A(a!) (5) 



with an arbitrary real function A(x), compensate for the change in the phase factor in equation (||), and the Wigner 
operator pw(P,X) remains unchanged. This ensures the gauge invariance of the definition (0). 

It is, however, quite clear that the choice of the integration path is still ambiguous and is not anyhow restricted 
by the condition of the gauge invariance. Another path of integration differing from the straight line would mean a 
differing definition of the variable P. This is the phase factor 



exp 



/ x+ y' 2 

ie 
ch 

\ X-y/2 



M*)**" (6) 



which would ensure the gauge invariance of the definition (||) for any integration path (Mandelstam, 1960). Vasak 
et al (1987) and Elze et al (1986) proposed removing the above ambiguity by choosing the straight integration path. 



1 Correct calculations using noninvariant quantities lead, nevertheless, to correct gauge invariant physical quantites (expec- 
tation values). The above mentioned references and section (v| contain examples of misinterpretations, approximations and 
assumptions that may result from attempts to provide physical reasoning to non-gauge-invariant equations, expressions and 
other terms used during the calculations. The discussion of the meaning of the slashed derivatives presented in subsection 
[II A implies that one should be cautious when considering the Poisson brackets for a non-gauge-invariant Boltzmann kinetic 
equation. 
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They showed that for this choice the transform (||) becomes equivalent to the transform (|l|) when substituting the 
covariant derivative d x + (e/c)A(X) for d . 

In spite of all that has been done, the situation remains rather confusing. The phase factor is introduced in equation 
(|J) only as a means of solving the problem of the gauge invariance. The arguments in favor of a specific phase factor 
are mainly intuitive. Hence, it is not quite clear which transformation is to be used if one prefers choosing other 
variables in the Wigner quasi-distribution function or if the dispersion relations for the charged particles differ from 
quadratic. It is our aim in this paper to derive systematically the transformation (|^) and to illustrate the importance 
and physical significance of the phase factor by several examples. 

Subsection II C formally defines the variables of the Weyl correspondence principle to be the kinetic momenta and 
the coordinates. We show that the form of the transformation (||) is due to the minimal coupling principle used to 
relate the canonical and kinetic momenta. The phase factor introduced in (Q) is shown to be a gauge transformation 
that transforms an arbitrary gauge into an axial gauge. Working with this axial gauge, the vector potential vanishes 
along the integration trajectory. In addition the kinetic and canonical momenta coincide and the two transforms, (jl]) 
and (H), are identical. We shall dwell on the connection between the kinetic and canonical momentum in the Wigner 
representation as well as on its classical limit. 

Quantum kinetic processes are in many cases treated by means of the Kubo linear response method Kubo (1957). 
An alternative approach, based on Keldysh (1965) diagrammatic technique, enables one to go rather easily beyond 
the linear response approximation. The first to formulate kinetic equations, based on the real time Green's function 
technique, introduced by Martin and Schwingcr (1959), were Kadanoff and Baym (1962). Fleurov and Kozlov (1978) 
using Keldysh formulation deduced two sets of linearized quantum kinetic equations which may in principle give an 
exhaustive description of the dynamic and kinetic properties of a system of interacting Fermi-particles. Since then a 
number of papers were published on quantum kinetic equations (Mahan, 1987; Rammer and Smith, 1986; Tugushev 
and Fleurov, 1983; Reizer and Sergeev, 1987; Khan et al, 1987; Reggiani et al, 1987; Davies and Wilkins, 1988; 
Bertoncini et al, 1989; Ferry and Grubin, 1995) and references therein. A gauge invariant formulation of the quantum 
kinetic equations for charged particles in an arbitrary classical electromagnetic field was derived in our paper (Levanda 



and Fleurov, 1994). Here, Section HI. briefly revisits this derivation and represents the resulting quantum Hamilton- 
Jacobi and Boltzmann equations in a more compact form by using the 'slashed' derivatives. These derivatives make 
the kinetic meaning of various terms in the equations more transparent. 

A hydrodynamic representation of the quantum theory was introduced by Madelung (1927) who transformed the 
Schrodinger equation, for a free particle, into particle and momentum continuity equations in the form similar to those 
found in hydrodynamic theories. We use our quantum kinetic equations in the same manner to obtain a hydrodynamic 
representation for description of dynamics of interacting charged particles in arbitrary external electromagnetic fields. 
The term hydrodynamic in this context may be somewhat misleading, as the reader may infer that 'hydrodynamic' 
approximations are implied. This, however, is not the case. The equations to be obtained would not assume any 
sort of approximations. The full set of equations for conditional moments (hydrodynamic equations) is equivalent 
to the original Dyson equations and, hence, they may account for all quantum effects. It is emphasized that the 
hydrodynamic equations in principle allow us to find the state of the system in and out of equilibrium without further 
assumptions (neglecting interactions these equations are equivalent to a single particle Schrodinger equation). This 
is not the case with other transport equations (e.g., for conventional quasi-classical Boltzmann equation). We show 
that collisions lead to an introduction of a 'collision force' in the momentum continuity equation. The collision force 
has the structure proposed by Bohm (1952) in his hidden variable interpretation, to modify the Schrodinger equation. 

Section [V] presents a few examples that emphasize different aspects of the Wigner quasi-distribution function and its 
equations. The Schrodinger equation and its solutions for an electron in harmonic potential and magnetic induction 
are found, in Subsection VA, from the hydrodynamic equation. This reveals the nature of the 'velocity' in such 
a system. We place a particular emphasis on the growth of the 'velocity' with the coordinates. The Wigner quasi- 
distribution function for such a system is found from the quantum kinetic equations. After integrating over the kinetic 
energy variable and taking one of the fields to vanish we reduce our result to expressions found in the literature for 
the Wigner quasi-distribution function in a magnetic induction or in a harmonic potential. The difference between 
the gauge-in varian t and non-gauge-invariant functions is discussed. 

Subsection |VB| formulates the Wigner quasi-distribution function for systems with discrete translational symmetries 
both in time and/or space (an example is an electron in a crystal). The acceleration theorem for an electron in a 
crystal is discussed and its relation to the gauge is pointed out. The Green's function of the heavy mass polaron in an 
electric field is found from the kinetic energy continuity equation. If the electric field vanishes, the Green's function 
is the same as the one found in the literature. If an electric field is present we find two solutions, one corresponding 
to a polaron with constant kinetic energy and the other one corresponding to a polaron with a diverging canonical 
energy. This serves as an example for the kinetic meaning of the slashed derivatives. 
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II. WIGNER QUASI-DISTRIBUTION FUNCTION FOR A CHARGED PARTICLE IN 

ELECTROMAGNETIC FIELDS 



An expression for a Wigner quasi-distribution function for a charged particle in an external electromagnetic field 
is derived. The coordinate transformation to the Wigner coordinates (^) is deduced rigorously from the Weyl corre- 
spondence principle. 



A. Derivation of the Wigner quasi-distribution function 



Weyl's correspondence principle suggests that an operator C, corresponding to a classical dynamical variable 
C(X, P can ) is obtained by replacing the variables P can and X in the Fourier transform 

C(y,vj) = — Lj J dP can dXe'^y p -+^C(P can ,X) (7) 

by the Hermitian operators, P ca n and X, representing these variables in quantum mechanics, 

C = J dyd^e^y p -+^C(y, m). (8) 

The expectation value of the operator (||) reads 

C) = J dx[p(j/,X)d] (9) 

where the angular brackets indicate an ensemble average. The second quantization allows one to present the density 
operator as 

p(y', X) = *t(x + y '/2)^(X - 2/72). (10) 
by means of the creation W(X) and annihilation ^(X) operators. Then its expectation value becomes 



^- <>/H)(yP ca „+™X) x 



C)= j dX' [^(X' + y'/2)CV(X' - y'/2)J> 
= J dP can dXC(P can , X) J dydme 

dX' (yi (x')eW h X yp ™ +ro * ) *(X')) (11) 
If the variable P can is chosen to be the canonical momentum, one can follow Moyal (1949) and write 

(c) = j^-^ J dP can dXC(P can ,X) J dydwe-^y p —+^ x 

/dX'(*t(X')e-»-/-e-/«e-..-/^ m ) 

= To *Tv4 / dP C a7i ( iX C (P can , X)py/can 

(Pcan,X). (12) 

where 

Pw(Pcan, X) = 

(pw(Pcan,X)) = J dye-^y p - (*t(X + y/2)*(X-y/2)) (13) 

is the Wigner quasi-distribution function. Properties and usage of it were discussed in the introduction. 

Equation (|b|) defines the Wigner quasi-distribution function, pw(P can ,X) which has, however, a serious drawback. 
Although, the physical quantity (C) is certainly gauge-invariant, the function p\y (P can , X), taken separately, depends 
on the gauge of the electromagnetic field. However, we are interested in a gauge invariant definition of the Wigner 
quasi-distribution function which can be derived by doing some modifications in the current derivation. 
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B. The need for a gauge-invariant Wigner quasi-distribution function 



A question which we find necessary to address is why a gauge-invarian Wigner quasi-distribution function is really 
needed? In the papers referred to in the Introduction a gauge invariant formulation was proposed only as a convenient 
representation, meant either to help the intuition or to directly correspond to a classical distribution function, with 
the coordinates and the kinetic momenta as variables in the classical limit of the Wigner quasi-distribution function. 
We want to emphasize that this is not the whole story a non-gauge-invariant Wigner quasi-distribution function 
does not always allow one to calculate expectation values of kinetic quantities. In particular we shall see that it is 
impossible to calculate the expectation value of the kinetic energy of a system in an electromagnetic field from a 
non-gauge-invariant function as phase-space integrals. 

For the Weyl correspondence principal to apply, the operator C must be a function of a linear combination of the 
operators X and P can (Moyal, 1949). Although, X and P can can be any complementary operators, in the following, 
we shall only consider the case where the two operators are the coordinate operator and the momentum operator 
(kinetic or canonical). The need for a linear combination of the operators can be seen by substituting C(aX + bP can ) 
for C(P can , X), with a, b being constants, in equation (Q), changing the variables to £ = aX + bP can , r\ = aX — bP can , 
7 = vj /2a + y /2b, fi — tu/2a — y/2b. 



((7 - r))bP ca „ + (7 + v)aX) 



d£dr; exp 



C(0- 



(14) 



Integrating over the new variables, rj, fi, 7 and then over £, one finds that C — C(aX + bP can ). The correspondence, 
C = C(X, P C an) goes into C = C(X , P ca n), which exists only for linear combinations of the operators. 

We conclude that the operators X and P ctm cannot be chosen arbitrarily and the Wigner quasi-distribution function 
should be set up for two non-commuting operators chosen according to the measured quantity (&^- I n general a 

kinetic operator, e.g., kinetic energy P 2 /2m, not always can be written as a linear combination of canonical momenta 
and coordinates, hence, its expectation value not necessarily can be calculated with a non-gauge-invariant Wigner 
quasi-distribution function. 

There are two cases in which equation ( |l2| ) can be employed. The first is that of a system in homogeneous and 
static electric and magnetic fields. Then the vector potential is a linear function of the coordinates. Hence, the kinetic 
momentum is a linear combination of the canonical momentum and the coordinate. The second case is when one is 
interested in the expectation value of the current (it is proportional to the momentum) . This can be expressed as the 
sum of the expectation values of the canonical momentum and of the vector potential, with each of these expectation 
values being calculated separately. 

The gauge-invariant Wigner quasi-distribution function will be formulated in terms of the kinetic momentum and 
will provide us a possibility of calculating expectation values of kinetic operators. 



C. Derivation of a gauge-invariant Wigner quasi-distribution function 

The Cartesian coordinates and the canonical momenta are not the only possible basis operators (Scully and Cohen, 



1987; Pimpale and Razavy, 1988). We may choose P to be the kinetic momentum then, using the equality (A.l 
derived in the Appendix A, equation (fill) becomes 



C = 



1 



{2nK) A 



dPdXC(P, X) 



1 



{2nh) 4 



dydrue 



-(i/H)(yP+vjX) 



dX' ( [&{X')#*x ) exp I ^ / yPA^X + sy)ds | *PO 



(2irh) 4 



where 



dPdXC(P,X) Pw (P,X). 



p W {P,X) - J dye 



(15) 
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exp 



/ , 1/2 

^ J y^(X + sy)ds ] (*t(X + y/2)*(X-y/2)) 

V -1/2 



(16) 



Here, as opposed to equation (|13|), the resulting Wigner quasi-distribution function, p\y(P, X), incorporates a phase 
factor and coincides with the definition (|J) postulated by Irving (1965). 

It is appropriate to discuss here the connection of the straight integration path with the axial gauge (see, e.g., 
Rohrlich, 1965). A vector potential A\(x) in any gauge can be transformed to the axial gauge by means of the 
transformation 



A x (x) -> A' x (x;x") 



A x {x)+dt 



I 



1/2 



V 

1/2 



\ 



{x"-xY J A^dx" + x)/2 + s(x" -x))ds 

-1/2 

1 



(x"-xY J dsq - s)F Xfl ((x" + x)/2 + s(x" - x)) 
(x"-xY 1 - e ' a j 1 + a) 



F x ^x") 



a— {x" — x) u B v 



(x" - xY Q - \{x" - xyP + • • •) F x ,(x"). 



(17) 



The resulting vector potential depends only on the electromagnetic tensor 

F Xfl (X) = dxA^X) - dpAxiX) 

and on the coordinate x" which can be chosen arbitrarily. This form of the vector potential has two important 
properties. In regions free of currents and charges, j X = (cp, j) generating the electromagnetic field, 



d X A' x (x) 



1/2 



{x"-xf J ds(--s) 2 j x ((x" + x)/2 + s(x" -x)) =0 

-1/2 



and 



(a:" 



X A' x (i 



0. 



(18) 



(19) 



Hence, this transformation converts any vector potential into a transverse vector field. 

Now substituting this axial vector potential in the definition ( [l6| ) and taking x" — x = a(x'—x)/2 = ay/2 with a being 
a constant, one sees that the integral in the phase factor becomes zero making the phase factor equal one. Therefore, 
the gauge-invariant Wigner quasi-distribution function coincides for this specific case with the gauge non-invariant 
definition (|l~3|). 



All that has been done thus far can be simply restated as follows. The choice of the axial gauge with x" 



ay + x, 



ensures that the vector potential indeed vanishes (see equation (|19|)) along the line that goes through x" and x. The 
displacement operator along this line, can equally be written as exp (— iP^ an (x" — x)^) or as exp (—iPkin( x " ~ x )^) 
and the Wigner quasi-distribution function of the canonical momentum variable coincides with the Wigner quasi- 
distribution function of the kinetic momentum variable. The function of the phase factor is to 'transform' an arbitrary 
gauge chosen for the vector potential into the axial gauge. 

Two remarks are in order. This derivation of the gauge-invariant Wigner quasi-distribution function, can be 
applied to other correspondence principles as well, which in turn will lead to other quasi-distribution functions. As 
was shown in references (McCoy, 1932; Cohen, 1966; Summerfield and Zweifel, 1969; Cohen, 1966, 1976, 1987), these 
quasi-distribution functions can be obtained from the Wigner quasi-distribution function by simple transformations. 

The equality derived in the Appendix A is applicable whenever the operators obey the relation 
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Pcan =PW+ A(X), 



(20) 



where Pw is the desired variable in the Wigner representation. That is, as long as one chooses a variable that obey 
(pp|), the transformation (|J) should be used. 

In the language of classical mechanics it can be stated that the relation (|2(]) holds if the coupling to the electro- 
magnetic field is via the term A(X)X in the Lagrangian L. Then the canonical momentum is of the form 

Pcan = 



dL(X,X,A) dL(X,X,A) 



dX 



dX 



del' 



A(X) = P W + A(X). 



(21) 



Whenever an effective Hamiltonian is used the relation ( po| ) is assumed to hold (see, e.g., Nenciu, 1991, and references 
therein). This however, is not necessarily the case. The particle may have a different dispersion relation due to its 
coupling to a set of degrees of freedom (e.g., polaron). Then if the Hamiltonian for such a composed quasi-particle is 
written, the minimal coupling (|2C|) may be not necessarily applicable. 



D. Kinetic and canonical momenta in the Wigner representation 

The phase factor in the Wigner quasi-distribution function, equation (|l6|), can be written as a Taylor series and 
integrated over s so that equation (Tlq) takes the form 



p w (P 1 X) = / dyexp 



My'dfWA.iX) 



tf{X + y/2)*(X-y/2) 



exp 



-Jo 



-d Pv d x 



A„(X)d p A x 



J dyexp (-^V p -) + y/ 2 )*( x - y/ 2 ) 

exp (- c An(X)d p ^ Pcan (P,X) 



(22) 



where $i.^(X) — jo(A/2)A u (X), A = Tid Pv d* , jo{x) — sinx/x al - x 2 /3 + • • •, and the derivatives with respect to 
X in the exponent act only on the vector potential. 

This is a general expression valid for an arbitrary electromagnetic field. We define pw(P,X) and P to be the 
Wigner quasi-distribution function and its momentum as derived from equation ( p2[ ) with the modified phase factor 
exp (f [A„(X) - A^X)] D p ») . The equality 



p w (P,X) = p w (P--A(X),X 
V c 



(23) 



then holds and one may say that the two momenta are related by the equation 



Pkin {X) = P(X) - -A(X). 



(24) 



If all the derivatives of the electromagnetic field of order higher than one vanish, meaning only static and homogeneous 
fields are applied, the relation (24) is simplified and P becomes the usual canonical momentum. The kinetic momentum 
and the canonical momentum as defined by the Wigner quasi-distribution functions in this section are related by the 
equality, 



can\x=x" 



Pkin(X) = P can (X;x") - (e/c)A(X;x") 



(25) 



where A(X; x") is in the axial gauge, (|T7|) . 
The differential version of equation (p5|) is, 



d*P kin (X) = d*P can {X;x") + dl P can {X;x"). 



(26) 
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Equation (|26|) shows us that the variation of the canonical momentum is composed of two parts: variation of the 
kinetic momentum due to kinetic effects and variation of the kinetic momentum due to the change of the gauge as the 
reference point x" moves. In particular we may use equation (jl^) and write (for systems in static and homogeneous 
field) 

^(X;x'0 = fE + f c ^xH=i W x"), (27) 

where (Lor{x") is the Lorentz force acting at the point x" . 

Let us imagine, a quantum system confined to a small part of the three dimensional space such that it can be 
described as being at rest (equivalently we may think of a heavy mass particle). Another classical particle with a 
well-defined trajectory is chosen to define the location of the point x" . The gauge part of the canonical momentum 
of the quantum system (as measured in the axial gauge by x" at the classical particle) will change in time according 
to the equation of motion of the classical particle (pfj). While all that time, the kinetic part of the kinetic momentum 
of the quantum system will not change since it is at rest. We shall return to this point twice, when discussing the 
heavy mass polaron and when discussing the 'acceleration theorem'. 

Now, suppose that the classical particle is another quantum system that interacts with the first system and exchanges 
momenta with it. Momentum conservation applies only to the canonical momentum. If the conservation laws are to 
be expressed in terms of kinetic momenta, the conservation laws should take into account the kinetic as well as the 



gauge part of the kinetic momenta. More on this point is discussed in Subsection V B 



E. The classical limit of the Wigner quasi-distribution function. 

The classical limit of the Wigner quasi-distribution function is discussed here. This limit has been already addressed 
by other authors (Heller, 1976, 1977; Berry, 1977) and here we give only the modifications due to the gauge invariance 
and appearance of an energy variable. A different approach was taken by Bund (1995). 

Following Van Vleck (1928) we express the quasi-classical limit of a wave function for a state a as 

rl> a (X) = N 1 ^ 2 V 1 / 2 e s/lh . (28) 

N is the normalization constant, S(X, a) = 5(R, t, a) is the classical generator of the canonical transformation from 
the P, R system to a set of new variables a, (3 and T> = det \d 2 S / dHda\ is the functional determinant. The Wigner 
quasi-distribution function reads 

p w (P,X\a) = 

d 4 y Vdet \d 2 S(X + y/2, a)/dRda\ det \d 2 S(X - y/2, a)/dRda\ x 

(29) 



exp 



- [(P - A{X)Yy„ + S(X - y/2, a) - S(X + y/2, a)] 



We neglected terms, non-singular in Ti, in the phase factor, since they vanish as fi — > 0. 

We conclude that the relation ( p5| ) is expected to hold in the quasi-classical limit as well as for quantum systems in 
static and homogeneous fields. In the classical limit the frequency and the wave length of the electron wave function 
are much smaller than any other time and space variations in the system, that is, the electromagnetic field is taken 
as static and homogeneous. This relation between the canonical and kinetic momenta are always valid in classical 
mechanics as well as for quantum mechanical operators. However, in the Wigner representation of the quantum 
mechanics this relation generally does not necessarily hold. 

III. DERIVATION OF KINETIC EQUATIONS. 

The diagrammatic technique proposed by Keldysh (1965) allows one to write the Dyson equation for a nonequi- 
librium quantum system in an arbitrary not quantized electromagnetic field. A semiclassical Boltzmann equation 
as well as its quantum versions may be obtained from this Dyson equation by making certain transformations and 
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approximations. It is our aim to present below quantum equations avoiding approximations as much as possible. In 
this sense, both the equations to be obtained and the Dyson equations are equally valid. 

Following the five step procedure presented in our paper (Levanda and Fleurov, 1994) we arrive at equations which 
are quantum analogs of the Hamilton- Jacobi and the Boltzmann kinetic equations 



2m 8m 



G(P,X) = 2a x 



exp 



a x S(P, X)G(P, X) + G(P, X)S(P, X)a x 



and 



G(P,X) 



(30) 



exp 



,_d 5_, 

dP dP 



V,£(P, X)G(P, X) - G(P, X)£(P, X)a Q 
where the following definitions are used 



- %eh (*) E(X) • V J 



= ^ + e i0 (f)E(X).V- 



f = P + ~eji (y) {- C B(X) XVP + E{XW 



t X =V x +ej 0[ - 



A ' ^B(X) x V p + E(A)d £ ) : 



(31) 



G — ( ] 5j — 



E a 



(32) 



j'o(a^) = sin(x)/cc and = sin(a;)/a; 2 — cos(a;)/a;. 

Here the derivative with respect to X in the operator A = hd p d x acts only on the electromagnetic fields E and B. 
The other derivatives, in the left hand sides, with respect to X and P act on all the functions. In the right hand side, 
the first and the second derivatives of in the exponent act on the first or on the second terms of the pair in the sum, 
respectively. 

This is a generalization of the so-called gradient expansion first used by Kadanoff and Baym (1962) for system 
without electromagnetic fields, by Tugushev and Fleurov (1983) for system in homogeneous and static electric field 
and by Mahan (1987) for homogeneous and static electric and magnetic fields. The form (^0|), ( |3l| ) is more general 
and valid for arbitrary electromagnetic fields. An alternative representation of the gradient expansion was given by 
these authors (Levanda and Fleurov, 1994) in terms of the phase loop function (see also Appendix B). 

These equations are equivalent to the Dyson equations for the Keldysh functions from which we have started. In 
fact, they are just alternative representations of the Dyson equations. Therefore they may provide an accurate and 
complete physical description for kinetic processes of a many particle system in an arbitrary electromagnetic field. We 
wrote two matrix equations which would correspond to six equations for complex functions. However, the number of 
really independent functions and, hence independent equations is only two (see, e.g., discussion in Fleurov and Kozlov, 
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1978). One can, for example, take an off-diagonal term in the matrix equation (|30|) in the representation (|32|) which 
would be an equation for, say, retarded Green's function G r determining the spectrum of the system. The second 
equation can be produced from one of the diagonal terms of the matrix equation ( |3l| ) (in the same representation) 
and this is the equation for the function G K which in the semi-classical approximation reduces to the conventional 
Boltzmann equation. 

Both equations can be expanded in powers of the derivatives of the electromagnetic fields (A expansion), without 
effecting their gauge invariance. This, however, does not correspond to the expansion in powers of h of equation (|3f]), 

The Liouville theorem in classical mechanics can be derived as a limit of equation (pll). To see it we neglect the 
interaction term (the right hand side) and integrate with respect to the variable e (the classical phase-space does not 
contain the energy variable) . Equation (pll) becomes 



d_ 

dT 



A 

ejo\ j 



E(X) + — P x -B(X) 

m c 



P V 



x 



G(P,X) = 0. 



(33) 



If all the derivatives of the fields of order higher than one, vanish, equation 
This should be compared with: 



3h becomes classical Liouville equation. 



1. The conditions for the coincidence of the Ehrenfest equations of motion for a wave packet in an electric field 
(see, e.g., Messiah, 1962, p. 220) with the classical equation of motion (if the electric field has at most a linear 
coordinate dependence). 

2. The conditions for a quantum system to be deterministic. Moyal (1949) showed that if the Hamiltonian of a 
quantum system is a second-degree polynomial in X and P (a uniformly accelerated particle or a harmonic 
oscillator), then the quantum system is deterministic in the sense that the time evolution of the Wigner quasi- 
distribution function is given by the classical Liouville theorem. That is, the quantum nature of the system 
enters only via the initial conditions (Bartlett and Moyal, 1949; Lee and Scully, 1983). 



A. Kinetic meaning of the slashed terms 



We shall now demonstrate the relation of the slashed energy, momentum and derivatives to the operators of the 
kinetic energy and kinetic momentum. This will clarify the physical significance of these quantities and their difference 
from the non-slashed quantities. We also give two examples of the Wigner quasi-distribution function that demonstrate 
the difference between the two derivatives. 

The slashed time derivative is defined as 



n G<{P,X) = 



d y exp I ——ry H — — 



te 
ch 



( 



d 4 y exp 



-rPy 



\ 



le 
~ch 



zftA + -A (X - y/2) 
ol c 



( 



d A y exp 



V 



1/2 

y^A li (X + sy)ds \ G<(y,X) 

1/2 

1/2 \ 

y"A fl (X + sy)ds 

1/2 / 

e 



-A (X + y/2) 



1 Py + 16 
h ch 



y^A ll {X + sy)ds ■ 



* ^ 6 A , X 

ih— + -Aq(x) 
at c 



c 

1/2 

/ 

-1/2 



c 
or c 



G < (y,X) 



(x')ip(x)) 



(34) 



This definition has been written for the G < (x,x') = ip[x,x') component. The expressions for the other components 
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are similar. The definition of the slashed spatial derivative is 



iht x G<(P,X) = I d 4 yexp - l -Py + -| j y»A^(X + sy)ds 



-1/2 

'(-ffiV 1 - ^A(x)) - (ihV x ' - ~A(x'))] i(ip + (x')ip{x)) 
In a similar way one can introduce definitions for the slashed energy and momentum 

1/2 



2iG<{P,X) = y"d 4 yex P _Ipy+i£ j y^A fi (X + sy)ds 



-1/2 



i(ip (x )ip(x)) 



(35) 



(36) 



2fG<{P,X) = / d 4 yexp 



/ 



V 



1/2 



h ^ ch 



y^A il {X + sy)ds 



-1/2 



-av 1 - -a 

c 



(37) 



The operators in the square brackets in (|36|) and in fl37j ) are the kinetic energy and kinetic momentum operators. The 
kinetic momentum operator that acts on the wave function of the system is gauge invariant and produces the values 
of the electron's mass times a quantity which has a meaning of an average velocity. 

We shall now show how the uniformity of the particle and the current density make y x G(P, X), but not 
V X G(P, X), vanish. The system is not assumed to be homogeneous, e.g., it may be placed in a non-uniform 
magnetic induction or in a time dependent electric field. That is, we assume only that the density and the current do 
not depend on the coordinate: 



tp*{x)ip(x) = p{x) = p 



and 



e 
2m 



4V 1 



A(xj) + (ihV x ' - -Afx'j)] iip*(x'U(x) = j(ar) 

/ V C /J x'^x 



and show that f x G{P,X) = 0. 

A single particle 'wave function' of the system is written in the form ip(x) = R(x) exp[iS(X)/K] and with the two 
aforementioned conditions for homogeneity, one finds, R(x) = yfp and 'V x S(x) — mj/e + A(a;). Putting this form of 
the wave function into equation (^5|) the derivative *f x Q{P, X) is readily seen to vanish. The derivative V X G(P, X) 
vanishes as well, only if the electromagnetic fields vanish. Substituting the same wave function in equation (|37j ) we find 
that J>G < (P, X) = mj/eG < (P,X) and PG<(P, X) ^ mj/eG < (P, X) if the electromagnetic fields have non-vanishing 
gradients. 

The Wigner quasi-distribution function for a charged particle in a static and homogeneous electric field in a system 
without interactions obey the equations 



h 2 



e - — + — ( V A + eEd 

Zm am 



Pw (P,X)=0 



(38) 



d_ 

df 



eEV' 



P 

rii 



(v x 



p w {P,X) = 



(39) 
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If the condition is assumed that the derivatives of Wigner quasi-distribution functionwith respect to the coordinates, 
V x p(P, X) = 0, vanish then the Wi gner quasi-distribution function reads 

PW (P,X) = 

{mheE)V3 Al ( ( TO ft e £)2/3 ( e ~ ^) ) ke=*? S ( £ ~ ^) ( 40 ) 

which is a time and coordinate independent solution, Ai is the Airy function. The solution ( f4(i| ) implies that the 
application of an electric field accelerates the particle and, hence, causes a broadening of the energy uncertainty. The 
broadening is of the order of (?iE) 2 / 3 and is only due to the acceleration, no interactions are assumed. The Wigner 
quasi-distribution function is a spectral function that yields an electron dispersion relation. It does not produce the 
dynamics of the system. The dependence on the energy and the momentum can be rewritten in the form 

P 2 /P eE ,\ (P-eET) 2 „ , N P 2 (T) 

e = e - eE — T T - ± '— = e - eER(T) v ' 



2m \m 2m J 2m 2m 

where R(T) and P(T) are the coordinate and the momentum of a classical particle that moves in the same electric 
field. This representation reveals the hidden acceleration of the particle. If a steady current flows in the system the 

slashed derivatives vanish and the dispersion relation found from equation ( |38| ) is S — , hence, no uncertainty 

broadening takes place. We conclude that the Wigner quasi-distribution function of a system without interactions 
that does not depend on the coordinates correspond to a particle accelerated in the electric field and not to a steady 
current system. 

Now turning to a system with interactions we may use at this stage the simplest mean free time approximation in 
the quantum kinetic equations. This approximation has been used for more than a century in classical and quasi- 
classical systems (see, e.g., Ziman, 1960). To introduce the mean free time we first consider a solution without any 
field p(P,X)\ E = 0. Then the Liouville equation ( p3| ) allows as to make the 'canonical transformation' 

e^e-eE- (R-Ro), 
P -> P - eET, 

P„ e 

to 



R -> R n —T + — ET 2 . 



m 2m 

in order to arrive to a solution in an electric field. The transformed Wigner quasi-distribution function after the 
canonical transformation has the following properties: (i) If the system, has been homogeneous prior to the application 
of the electric field, it remains such afterwards; (ii) The slashed derivative vanishes. 

The next step is to assume that after the interaction converts the variable T into a parameter having the meaning 
of the mean free time that needs to be determined from the quantum kinetic Boltzmann equation. The mean free 
time may be a function of the momentum and the energy, r(e, P). The derivative with respect to T in the kinetic 
equations vanishes (T is now a parameter) as well as the derivatives with respect to X and e (due to the canonical 
transformation of e). In the left hand side of the Boltzmann kinetic equation ( |39] ) we are left with only one term, 
eE • V P /9u/(e, P, R), just like in a quasi-classical Boltzmann equation (see, e.g. Ziman 1960). This term should be 
equal to the collision integral in the right hand side of the equation, the electric force is equal to the collision force. 
In the left hand side of the Hamilton-Jacobi kinetic equation ( |38|) we are left with [e — P 2 /2m)pw(s, P, R) that is 
equal to the interaction which causes a broadening in the dispersion relation and vanishes for classical systems. 

Let us look at the dependence of the Wigner quasi-distribution function on the mean free time. The Wigner quasi- 
distribution function in equilibrium is usually written in the form pvK (P, X) | E=0 = A(P, X^npis) where np(e) is 
the Fermi function and A(P, X) is the spectral function. In the linear approximation the collision integral does not 
vanish only if nonequilibrium corrections to the Fermi function are considered, the spectral function may be taken 
at equilibrium (Fleurov and Kozlov, 1978). Hence, the linear correction to the Wigner quasi-distribution function is 
taken in the form 



PW {P, X) = A(P, X) [n F {e ~ eE • (R - R ))] 



n F (e)- f^r) eE-(R-Ro) 
= pw(P,X)\ e=0 + A(P,X) (^fj ^± T ( £ , P ,R). (41) 

This form of the Wigner quasi-distribution function may be inserted into the quantum kinetic equation in order to 
calculate the mean free time r(e,P,R). 
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IV. KINETIC MOMENTUM CONTINUITY EQUATION 



Quantum kinetic equations provide us a possibility to calculate Green's functions. Green's functions, however, are 
not directly measurable, and are only used at certain steps of calculations. Green's functions enable one to calculate 
quasi-distribution functions, and macroscopic measurable quantities are calculated as their averages over the phase 
space. Although deduction of various equations for the macroscopic measurable quantities from the quantum kinetic 
equations is possible, we shall consider here only the kinetic momentum continuity equation. The particle continuity 
equation and the quantum mechanical equivalent of the classical Hamilton- Jacobi equation were discussed in Levanda 
and Fleurov (1994). 



A. Space conditional moments 



Macroscopic measurable quantities can be written in terms of gauge invariant space conditional moments, P n fJi (X) 
which can be defined as the averages of various powers {P^) n of the momentum at a given point, X, 



1 



1 



p{X) {2ttK) a 
d 



d 4 P(P^ n p w (P,X) 



-ih- 



dy 



1/2 



Here 



exp 



P(X) 



ch 



y»A^X + sy)ds p(X,y) 



-1/2 



(42) 



d 4 Pp w (P,X) 



(2irh) 4 

is the particle density. The canonical (non gauge-invariant) space conditional moments are defined as 

P 71 can,/^(^0 



11/* 1 f \ n 

p(x)(2W J ^WW™*) = pjx) {- lh o-y) piX ' 



v) 



(43) 



The relation between the canonical space conditional moments and gauge-invariant space conditional moments is 
trivial if equation (|25|) holds. In the general case one can also see that the relation 



P(X;x") = Pcan(X) - {e/c)A(X;x") 



(44) 



holds regardless of the nature of the electromagnetic field (as well as equation (EJ)). 

This can be readily understood, since P C an(X) is merely the expectation value of the current operator times the 
electron mass divided by the density, 



-ifi 



Pcan{X) 

[d x ' - d x 



_ mJ can (X) 
ep(X) 



^(X')^(X) 



(45) 



In a similar manner we have an expression for the gauge-invariant terms 



P(X) 



mJ(X) 

ep(X) 



1 (-ihd x ' - lA{X') + ihd x + fA(A)) ^(X')^(X) 



*t(X')#P0 



(46) 



X'^X 



13 



Relations between the canonical and gauge- invariant conditional moments are best seen with the help of the cummulant 
expansion generated by the function 



log 



/ 1/2 



exp 



- J y^A fi (X + sy)ds | p(X,y)/p(X) 

V -1/2 



1/2 

= ^ J y^A^X + S y)d S + \og(p(X,y)/p(X)). 

-1/2 



Then writing the density function in the form 



p(X, y) = p(X - y/2)p(X + y/2) exp — (S(X - y/2) - S(X + y/2)) 



one obtains cummulants. 



and 



«2n+l 



2n+l 



1/2 

| / y»A M (X + 8y)ds 



-1/2 

-ih d \ 2n 
2 dX) 



■Ax) 



«2n — 



\ In / n \ 2n 

-in\ I a 



\dX 



io g p(x)+ -in 



'dy 



2n 



1/2 



ie 
ch 



y^A^X + sy)ds 



-1/2 



/ \ 2n / o \ 2n 

-in\ I a 



\dX 



logp(X). 



n = 0, 1, 2 . . . The odd cummulants contain canonical cummulants plus cummulants of the phase factor. The asymme- 
try of the Wigner quasi-distribution function is, hence, represented by these odd cummulants, therefore the asymmetry 
of the gauge-invariant and canonical momenta are different. The even gauge- invariant and canonical cummulants are 

the same and, in particular, the conditional mean square deviation of the momenta, K2 = P 2 — P , is the same in 
both cases and the uncertainty relation holds for both the canonical and the gauge-invariant momenta. 



B. Derivation of kinetic momentum continuity equations 

Here we discuss an equation which can be called a kinetic momentum continuity equation. It is obtained from the 
trace of equation ( |3l| ) , after multiplying it by P and integrating over the energy and momentum, 



PxP(X) 



in 



-eE 



(X)p(X) - (e/mc) (p(X)P(X)) x B(X)- 
^/dPP W P,X). 



(47) 
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The kinetic momentum continuity equation ( pl7j ) contains a collision integral, multiplied by P and integrated over 
the momentum space. It describes the transfer of electron kinetic momentum from and to a four dimensional volume 
element, at X, caused by the collision processes. 

Equation ( |47| ) can be rewritten in a more transparent form 



where 



dv(X) 
Tfi — — : — 
dT 



e— , 



eEiX) - -W(X) x B(I) - £v^^ + f co „(X) (48) 



fcoii(X) = 1 | dP (P - mv(X))I coll (P,X), 



v(X) = P(X)/m = 3 km {X)/ep(X) 



and the equality 



(P X P(X)-P X (X)P(X) 



(49) 



is used. The left hand side of equation Q4SQ is the gradient of the so-called 'quantum potential', whereas the term 
o~\v = {h 2 /4m)p(X) V A V In (p(X)) is the stress tensor used in the hydrodynamic formulation of the Schrodinger 
equation (Madelung, 1927; Holland, 1993). The electron density continuity equation (equation (47) in Levanda and 

Fleurov, 1994) shows that the quantity ep(X)P(X)/m is the kinetic particle current density Jkin(X) at the point 
X, hence, v(X) is a quantity having a meaning of an average velocity of the particles. This averaging includes both 
thermodynamical and quantum averagings. 

The collision term f co u(X) appears due to possible processes of electron scattering. Bohm (1952) (his equation (31)) 
argued that hydrodynamic equation in the quantum theory should contain terms depending on the difference p — rav 
but did not present a physical mechanism that may lead to such a term. It is clear that the collision term f co j;(X) is a 
function just of this type. Its specific form is determined by the collision integral and can be now calculated explicitly 
for any particular system. 

This continuity equation has been obtained from the gauge invariant quantum kinetic equations which depend 
explicitly only on electromagnetic fields and do not depend on their potentials. These equations are actually Dyson 
equations for the Keldysh Green's functions in a differential form. In the absence of collision processes, they can be 
reduced to a Schrodinger equation and to a quantum generalization of the Liouville theorem (Moyal, 1949). All the 
physical processes described by the Schrodinger equation can be described by the quantum kinetic equations without 
any need of using potentials or vector-potentials of the field. The continuity equation for the velocity v(X) can be 
looked at as a quantum generalization of the second law of Newton. Equation (48|) converts into a classical Newton 
equation of motion with a Lorentz force, if the quantity v(X) is replaced by the classical kinetic velocity and h — ► 0. 

We mention, here, a few general things about the collision terms that emerge from their explicit expressions derived 
from equation (38) of Levanda and Fleurov (1994). The collision force is written in the form 



p(X)f coll (X) 



-ih-V y D(X, y) 



mv(X) + -A(X) D(X,y) 
c 



(50) 



where 



D(X,y) 



yo=o 



d 4 Xl 9 H (X -( Xl ) ) D{X,y\ Xl ) 



yo=o 



and 
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b{X,y\ Xl ) = [S>(X ,X-y,( 2 ; 1 ) ,x 1 )G < ((x 1 ) ,x 1 ,X ,X + y) 

Vo=0 

-S < (X ,X - y, (a;i)o,Xi)G > ((xi)o,xi,Xo,X + y) 
-G > (X , X - y, (xi) , Xi)S < ((xi) , Xi, X , X + y) 
+ G < (X ,X-y,(x 1 ) ,x 1 )£ > ((x 1 ) ,x 1 ,X ,X + y)] . 

If the Green's functions and the mass operators obey the equalities G K (X, Xi) = — G > (a;i,X) and S < (X, xx) — 
— E > (xi,X), then D(X, y\xi) = D(X,—y\xi) and all the processes in the system are detailed balanced as discussed 
in Section 4 of Levanda and Fleurov (1994). In particular creation and annihilation of particles, D(X,y) = 0, 

and the collision force, V y D(X,y) = 0, vanish. 

y->0 

If transport processes occur in the system, the detailed balance no longer exists and the conservation of the number 

of particles or a vanishing collision force are not necessarily true. The function D{X, y\x\) = can be regarded 

as describing the transfer of particles to a four dimensional volume element, at X , from a four dimensional volume 
element, at xi, caused by the collision processes. If the net transfer of particles to a point due to collision processes 
does not vanish, the momentum at that point changes by the amount equal to the mass transferred in a unit of time 
times the longitudinal part the velocity (the velocity without the rotation part). This is the meaning of the term on 
the right hand side of equation ( |50| ) . 

In the same way h'V y D(X, y\x\) stands for the transfer of momentum from a four dimensional volume clement, 

at X, to a four dimensional volume element, at x±, caused by the collision processes. It is the quantum mechanical 
equivalent of the force applied to part of the system. This is the meaning of the term in the right hand side of equation 

(p0|). It obeys Newton third law, TiV" D(X,y\x\) = —K^ y D(xi,y\X) only for internal interactions. If for 

. . y ^.°. . . y ~*°. 

example an interaction with a phonon bath or impurities is considered this relation will not hold. 

It has long been known that approximations used in transport theory can in principle lead to transport equations 
that are not necessarily consistent with all the macroscopic conservation laws (see, e.g., Mermin, 1970; Hu et al, 
1989). This possible inconsistency was sometimes used to estimate the correctness of the results obtained (Thornber 
and Feynman, 1970). It is therefore desired to build the macroscopic conservation laws into the structure of the 
approximation used. Papers by Baym and Kadanoff (1961) and Baym G and Kadanoff (1962) (see also Green et al, 
1985) found criteria for maintenance of the macroscopic conservation laws in systems with electron-electron interaction. 
The main motivation for including the conservation laws in the transport formalism comes from the division of the 
transport processes into two types. First, the particles interact via the interaction specified in the approximation, with 
this process generally having a small time constant. Then, the kinetic process proceeds, with a large time constant, 
via interactions with an unspecified bath. This part of the interaction is described by the macroscopic conservation 
laws. Using the hydrodynamic equations it is possible to build the microscopic conservation laws into the structure 
of the approximation. That is, equation ( f48| ) holds for every X. This is a new aspect of transport theory that comes 
out naturally from the hydrodynamic equations. 

As an example, we mention the calculations of the conductivity in mesoscopic structures based on the quasi-classical 
Boltzmann equation. In these calculations the integration over the wave functions always precedes the evaluation of 
physical quantities (e.g., current density). Only after the integration, one equates the total influence of the electric 
field to the total influence of the impurities or of the phonons or of any other interaction. This procedure leads to 
macroscopic momentum conservation in the system but the microscopic momentum conservation is abandoned. In a 
paper to be soon published, we use the hydrodynamic equations to look for the current density in mesoscopic structures 
keeping the microscopic momentum conservation. We found that for elastic impurity s-scattering the current is close 
to being constant over the cross-section of the structure (i.e., the velocity of the electrons is higher near the edges of 
the structure) . This result differs from the commonly used results that assume constant velocity. We deliberately used 
the simplest and most widely applied approximations that are well known and commonly used for decades with the 
quasi-classical Boltzmann equation to emphasize the role of the one particular assumption we changed. We require a 
microscopic balance of forces rather than a macroscopic balance of forces. 



C. A single particle Schrodinger equation for a many particle system 

The hydrodynamic equations were originally derived from the Schrodinger equation without electromagnetic fields 
by Madelung (1927) and with electromagnetic fields by Janossy and Ziegler-Naray (1964). Here we discuss an interest- 
ing possibility. One may introduce a modified Schrodinger equation for a wave function of a single particle excitation 
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in a multi-particle system which is equivalent to the above hydrodynamic equations. We mean that it is possible 
to derive hydrodynamic equations for a system in electromagnetic fields with interaction, (p!8|), from the following 
modified Schrodinger equation, 

dT 



T n 2 



mv-(e/c)A(x) + r dT% oll (x,r) 



2m 



eip(X)- 



| X dX'f < UX , ,T) + J d P/ coH (P,X)^ 9{X). (51) 

Here i* oll (X) and f l coll (X) are the transverse and the longitudinal components of the collision term i co u(X). 



Now we demonstrate the derivation of equation (Mq) from equation (51) in two possible ways 



1. Introducing p(X) = and the Hamilton principal function S(X) = (-ih/2)ln(^(X)/^*(X)) (that 

is writing ^(X) = y/ p(x) exp(— iS(X)/K)) one finds the continuity equation and the expression for the canonical 
momentum 

¥ can (X) = V • S(X) = 
(mV(X) + (e/c)A(X) + ^ dT'f* oH (X, T'^j . (52) 

It should be noted that the velocity is derived from the wave function by means of the equation 

mJ kin (X) _ 



p(X)v(X) = 



2 I c 



A(X')- | T dT'f* o; /(X',T') 



^-^)-/ r ^4 w (x,T0)(*t ( x0^(x) 



(53) 

X'^X 



In the absence of the interaction terms this equation converges to the commonly used expression. 

Now differentiating equation ( |52"| ) with respect to time and using ( |5l| ) we arrive at the kinetic momentum 
continuity equation (f!8|). 

2. The same Schrodinger equation can be derived in a different manner. The expectation value of the kinetic 
energy density 



e{X) = -^rj d*PE Pw (P,X) 



P{X) 

satisfies the quantum mechanical equivalent of the classical Hamilton-Jacobi equation 



mv(I) 2 h 2 ( V x p(X)Y 



8m \ p(X) 



£int{X) 
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mv(A) 2 , h 2 V X2 ^pJX) 

'Si n t{X). (54) 



2m 



The last equality is obtained by integrating by parts in order to show explicitly the role of the quantum potential 
in changing the kinetic energy. This step is legitimate since an observable is obtained after integrating over the 
spatial coordinates. Substituting ( |52| ) into equation d54| ) the expectation value of the total energy of the system 
is expressed as a functional 



(E) = J dXp(X) (e(X) + <p(X)) 



J dX**(X) 



( - (e/c)A(X) + f dTft oll (X, T') N * 



V 



2m 



+ Lp(X)+£ int (X)\V(X), (55) 



where <p(X) is the electric potential. The Schrodinger equation is obtained by varying equation ( p5| ) with respect 
to the wave function ^>(X). This equation is, in fact, written for the wave function of a one particle excitation 
in a many particle system. It contains a complex potential £i n t(X) created by all possible scattering processes. 

Its imaginary part is simply Ime~i nt (X) — (h/p(X)) J dPI co ii(P,X). These sort of terms are discussed, e.g., 

within the framework of the optical model (see, e.g., Hodgson, 1963). The real part, Re£,- n t(Jf) is equal to 
/ dX f l coll (X) plus linked cluster expansion terms. 



D. Rotation of the velocity 



The phase of a single particle wave function represented by the Hamilton principle function S(X), must be single 
valued and satisfy certain boundary conditions. In a collision free system the boundary conditions lead to the 
Aharonov-Bohm effect (Aharonov and Bohm, 1959) and their controversial physical role has been discussed previously 
in several papers (Bohm and Hiley, 1979; Wodkiewicz, 1984; Liang and Ding, 1987; Peshkin and Tonomura, 1989). 
It has been argued that the hydrodynamic formulation uses field strength, rather than potentials, and, hence, it can 
describe the Aharonov-Bohm effect only when boundary conditions are added (e.g., Philippidis et al, 1981). Here, it 
will be shown that the Aharonov-Bohm effect can be derived from the hydrodynamic equations themselves without a 
use of potentials. The effect will be generalized to account for the role of collisions and the boundary conditions for 
the velocity in systems with interactions will be presented. 

To find the boundary conditions for the velocity v(A), its rotation is calculated. The rotation of equation ( |48| ) and 
the identity 



dB <9B 

- = - + Vx(Bxv) + v(V.B) 

yield the equation (compare to the derivation of Lorentz force from the Faraday's law, e.g., by Jackson (1975) 



(56) 



-eE(X) - -v(A) x Bpf) + f coU (X) 
c 



dl 



^ J B(A)-ds + ^f coH (A)-dl. 



(57) 



Equation ( p7| ) can be represented in the form 

m l V(X) • dl - m I V(X) ■ dl 



B=f co!! =0 



dT i> l coll {X)-d\ 



J B(X) ■ ds 



(58) 



where the integration in the left hand side of the equation is carried out along the curve bounding the integration 
surface in the right hand side. 
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The term § (v(X)) -dl| B=f _ Q is the rotation of the velocity for the colhsionless system without a magnetic 

induction. It can be calculated by means of the continuity equation and equation ( fl8|) which for a steady (dv/dT = 0), 
cylindrically symmetric system, coincide with the Schrodinger equation for particles with a well defined angular 
momentum. The Bohr-Sommerfeld quantization rule, therefore, follows 



jmv{X) ■ dl 



= nh 

B=f coi ,=0 



with n being an integer. A change by one quantum level in the rotation of the velocity corresponds to a change of a 
flux unit $o — ch/e — cl^hje. 

Equation ([57]) can be written with the help of equation ( |52| ) in the form 

~ j> V • S(X) = 2nn, (59) 

which is, in fact, the condition that the wave function is single valued. 
Equation ( |58| ) is rewritten in the form 

m i^{X) -dl = nh+~J B(X) ■ ds + J dT <j> i coU {X) ■ dl. (60) 

It states that the sum of the angular momentum of the electron, of the magnetic induction and of the system, that 
applies the collision force, should be quantized. All these three sub-systems are coupled together and a continuous 
change in the angular momentum of one of them results in a continuous change in the angular momentum of the 
others. For example, if the magnetic induction is raised gradually to its finite static value, it induces an electric field 
which changes the velocity lines in such a way as to keep the total angular momentum constant (Lentz law) (c.f., 
Peshkin, 1981). Equation (^0|) tells us nothing about sudden changes of the total angular momentum in quanta of 
h. This is done via coupling to an external system that is not included in our description (e.g., a current source that 
drives the solenoid). The external system that applies a torque to change the angular momentum of the electrons 
and the magnetic induction can change their angular momentum only in quanta of h, hence, if it applies a small and 
uneven torque to the magnetic induction and to the electron, it can change the angular momentum of the electron 
but not the total angular momentum. The conserved angular momentum of the external system in this interaction 
does not mean that this is a forceless interaction. The internal distribution of angular momentum in the external 
system has changed. 

Although the change of the velocity of the electrons is due to the induced electric field, the system does not relax to 
its state before the fields were applied just as a classical particle that does not relax to its initial velocity after the field 
has been removed. Relaxation may occur only if a dissipation mechanism is introduced into the system, via f co u(X). 
Equation (|^) explicitly shows that the change in v(X) can be viewed at as being due to a local action of the electric 
field at earlier times, or as a non-local action of the magnetic field at the present time. The collision force can cause 
a non-zero rotation of the velocity in the same manner as it can be done by an electric field generated by a time 
dependent magnetic field. 



V. EXAMPLES 



A. Wigner quasi-distribution function for a harmonic potential in a permanent magnetic field 

1. Schrodinger equation 

A system with a harmonic potential mw 2 z 2 /2 in the z direction and a magnetic induction B = (cm/ e)uj h in the 
x, y plane is considered. Before the application of the magnetic induction the electron is moving with the momentum 
p° = rav in the x, y plane. After its application the electron velocity changes according to equation (^8|) and becomes 

v(z) = p°/m + (e/mc)B x z = p /m + ujh x z. (61) 

The force acting on the electron stems from three sources, harmonic potential, magnetic induction and the quantum 
potential. Putting them all in equation (^8|), and assuming that the system is in a steady state, we find an equation 
for the balance of forces in the z direction, 
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/ 2 2 x h 2 d I d 2 y/p{z) 

2m dz \y/p(z) oz 2 



ill I ( l~ i / /II " \ 

(62) 



After integration (here, a particular choice of the gauge is made taken, see equation ^J) we obtain a time independent 
Schrodinger equation 



fi 2 d 2 ^p n (z) m{u 2 c +u> 2 H )z 2 

°-2^^^ 2 P X " HZ + E " 

h 2 d 2 ^p n (z) m(u: 2 c +uj 2 H ) , 2 2 



2m dz 2 

Here 



{(z + z ) 2 - (z a ) 2 ) + E n . (63) 



p°xu)jj 

z — 



m{u 2 H +u3 2 c y 



is the deflection of the center of the effective harmonic potential due to the Hall effect (it is the point where the Hall 
force equals the harmonic force). E n and p n {z) are the eigenenergy and the electron density for the nth eigenstate. 
The Hamiltonian can be written in the form 



P 2 (p° x ^ EL 



2 



. 2 

£T= ±- ?J- + —l— + —^-[2 + p<>x 

2m 2 



p 2 — I p° X 



2 



v^LL + m^ H + ^) {2 + z0)2 (64) 



2m 2 

If only the harmonic potential is present or if only the magnetic induction is present, the above Schrodinger equation 
takes a form similar to that which can be found in many text-books (see, e.g., Landau and Lifshitz, 1977). 
Solutions of this equation read 

VPn(z) = 



11, \ JjJlj*+u 2 h (z + z°) 



where H n are the Hermite polynomials. The energy levels are given by 

(p° x u H y 



E = {n + l/2)hduJ 2 c +LJ 2 H 



m{uj 2 c + u 2 H ) 



and, as expected, correspond to the levels of the effective harmonic potential shifted due to the Hall force. 

The above derivation emphasizes a connection between the momentum continuity and the Schrodinger equations. 
An attention should be also drawn to the velocity profile. It is the velocity (but not the current) which grows with 
the coordinate z. This can be seen directly from equation (^l|) or (more roughly) from the independence of p(z) from 
the coordinates x and y. This suggests that the canonical momentum is constant, whereas the kinetic momentum 
equals to the vector potential that generates the magnetic induction. 
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2. Wigner quasi- distribution function 



Gauge-invariant Wigner quasi-distribution function should satisfy the quantum Hamilton-Jacobi and Boltzmann 
equations for the system under consideration. The quantum Hamilton- Jacobi equation is found by substituting the 
electric field E = mu 2 z and the magnetic induction B = (cm / e)uj h into equation (BT 



1 

2irh 



;2 r 



muj 2 zd £ + tiiijJh x V p 



1 

m 



z h 2 

p- -— mutdP'd* 



pw(s, P,x) = 0. 



(65) 



Quantum Boltzmann equation (pl|) has the form 



1 

2m 



z 12 



V x - mu 2 zd e 



mu)H 



x V p 



1 

2m 



V - muj 2 zd e + mu> H X V p 



z?i z 
z 12' 



0. 



The set of functions that solve these equations are 




2m j 

where L n are the Laguerre polynomials. Q is defined by means of the equation 



1 



1 



= n 2 +e-\-+n) h\/u* + u 2 + -mujtz 2 - ^- 

2 2m 



iPx ~ muj H z) lo h 



2m uj 2 + u) 2 



H 



and 



Z n = 



Integrating equation (|67|) over e we find 



p w (P, z\n) = exp 



(p z - muj H z)uj H 
m {u) 2 c + u3 2 H ) 



m{ul+u 2 H ) 



'H 



L r> 



m(u> 2 



(66) 



(67) 



(68) 



Equation ( |6q ) is now compared with other existing expressions for the Wigner quasi-distribution function appearing 
in the literature for a harmonic oscillator (loh — 0) or for a particle in a magnetic induction (u) c — 0). 

The Wigner quasi-distribution function and an equation for it in one dimensional harmonic oscillator has been 
derived many times in the past (e.g., Heller, 1976; Bartlett and Moyal, 1949; Groenewald, 1946; Takabayashi, 1954; 
Carruthers and Zachariasen, 1983). All these derivations were done for a non-gauge- invariant Wigner function and 
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without the energy variable. The commonly used equation for the non-gauge-invariant Wigner quasi-distribution 
function for a harmonic oscillator is 



8m 



d 2 2 2 d 2 



= 0. 



p2 

2m 



— —i mu z )■ pw(P, z \ n ) 



(69) 



This equation does not converge to equation fl67| ) integrated over e. However, the solution of ( p8[) satisfies both 
equation ( |69| ) and (65) integrated over e. This situation of one solution for two different equations, becomes clear if 
we remember that e can — hu>(n + 1/2) and that 



J de ep w {e,P,z) = 



7 * P 2 1 2 2 

— aicj 1 — muj z 

12 6m 2 



J dep w (e,P,z). 



We conclude that although the Wigner quasi-distribution function, after integrating over e, is identical in the 
gauge-invariant and in the non-gauge-invariant cases, the equations they obey are different. This is due to the 
different expressions for the kinetic and canonical momenta. Before the integration over e, the two functions are 
different as their equations are. The fact that the two solutions are identical is connected with the fact that the gauge 
used in the non-gauge-invariant formalism assumes a vanishing vector potential, hence, integrating over e has washed 
out all the traces of the gauge. 

The Wigner quasi-distribution function for an ensemble of oscillators at temperature fc^T —1/(3 can be found using 
the linearity of the distribution function. 



p w (e,P,z\ > 



E 



3 -(n+|)ftW/9 



p w (e,P,z\n) = 



-(n+i)fiw/3 



1 m^ + wj. 



exp 



ly/uJ 2 c +U)' 2 H 



irk 



{z + Z ) 2 + 



itTi 
6ft 2 



in 2 (p /24ft 2 

6mw 2 ft 2 I h 



raujt. 



L, 



m(u> 2 c + u> 2 H ) 2 



iz + Z ) s 



Pz 

2m 



(70) 



Integrating over e, using the equality (see, e.g., Gradshteyn and Ryzhik, 1980, equation 8.975) 

1 



1 



^r = ]TL n (&K, |o|<i 



and normalizing, we arrive at an expression that converges to a well known expression (see, e.g., Hillery et al, 1984) 
for harmonic oscillators (at u>h = 0). 



pw(s, P, z|/3) = tanh 



h^/u> 2 c + U3 2 H (3 



exp 



UJ 2 + u> 2 



: tanh 



Tisju 2 + u 2 H (3\ fm{u3 2 c +u 2 H 



H 



-(z + Z Q ) 2 + V 



,2 

2m 



A comparison for the zero magnetic induction case can be done after substituting uo c — in equation (68). 
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Pw(P, z\n) = exp 



( 4 pI+pI 



exp 



2 P 2 x +pI 



Tiujh 2m 



" \hu)ff 2m 



(71) 



where 



2 +n j^-aS- 



This is the expression used in the literature (see, e.g., book of Mahan ,1987 | 



B. Wigner quasi-distribution function in a periodic system 

It is well known that the concept of a quasi-momentum plays an important part in the theory of electrons in periodic 
potentials (see, e.g., Weinreich, 1965). If the potential is also periodic in time, then a quasi-energy is to be introduced, 
Zel'dovich, 1973; Zel'dovich et al, 1976). Hence, if a periodic (in space and/or in time) electromagnetic field is applied 
to otherwise homogeneous system, it may be reasonable to use a Wigner quasi-distribution function formulated in 
terms of quasi-momentum and/or quasi-energy. This may provide us a better physical insight and make calculations 
easier. Here we shall try to reformulate the Wigner quasi-distribution function using the so-called fcg-representation, 
introduced by Zak (1972), assuming periodicity in space and/or time of the electromagnetic field. This may shine a 
new light on the 'acceleration theorem' related, as we shall see, to gauging rather than to electron dynamics. 

A /cg-representation can be constructed for any pair of operators A and B that satisfy the commutation relation 
[A, B] = —ih. Then a four dimensional constant vector a M , which is the space-time period of our system, allows one 
to define the operators 



T(a) = exp(iB A 'a Al ) and f(2ich/a) = cx_p(iA* 1 2Trh/a fJ _). 
Eigen-functions of these operators satisfy the conditions 

f(a)<p kq (x) = exp(ikv afl )(p kq (x) 



(72) 



(73) 



f(2-Kh/a)ipk q (x) = exp(iq ,1 (2'!T?i/a) ll )ipkq(x). 



(74) 



and form a basis of the fcg-representation. 

The operator A is chosen to be the coordinate operator x — (t,r), and B to be the energy-momentum operator 
p = (e, p) . Then k is the quasi-energy- momentum and q is the quasi- four-coordinate. Both the coordinates and 
momenta, as well as their corresponding operators, can be written as (Zak, 1972), 



and 



x = q + (m^a^a = q + X 
p = k + (m^(2Tih/a) lJ ){2TTh/a) = k + K 

x = q + X 
p = k + K. 



(75) 



(76) 



The 'correct' and 'transport' Green's functions used by Mahan are actually non-gauge- invariant and gauge- invariant forms 
of the Green's functions. 
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Here is a vector of integers, K is the momentum conjugate to the quasi-coordinates and X is the coordinate 
conjugate to the quasi-energy- momenta. It is clear that at h/a — ► the spacing between the discrete values of K 
vanishes and one approaches the continuous limit, K — > p, q — > x. The opposite limit h/a — > oo also yields k — > p and 
.V • x. 

The eigenfunction of the operators T{2izh/d) and T(a) can be now represented in the form 



(2tt)4 



2J exp (ik> J 'a fi (m x a\)) <5 (x M - g M - (m A aA)a M ) , (77) 



where Vo is the volume of the four dimensional periodic unit cell. These eigenfunctions can be used to define the 
creation and annihilation operators ofi (k, q) , a(k, q) in the fcq-representation in terms of the creation and annihilation 
operators tp>(x), ip(x) in the x representation, 

<p(x) = J dkdqa(k,q)ip kq (x). (78) 

a{k, q)=J2 e- ik " a ^ mXa ^if (q, + (m x a x ) atl ) (79) 

Transforming the Wigner quasi-distribution function into the fcq-representation the matrix elements 

(kq \p\ k'q') = -ih(d/dq)S{q' - q)5{k' - k) (80) 

and 

(kq \x\ k'q') = \ih(d/dk) + q] 5(q' - q)S(k' - k). (81) 

are used. Then the equality 

(*t(x + y/2)*(X-y/2)} = 

J dkdq(aHk,q)expUy^dUa(k,q)y (82) 
is obtained and equation (|l3|) in the fcg-representation now reads 

P w(P,X) = [ dye-^ h ^ f dkdq(J{k,q + y/2) a {k,q-y/2)) , (83) 



as can be verified by substituting equation (|79[). 

One thing that can be learned from equation ( [83] ) is that if all the particles in the system acquire an additional 
quasi- momentum 5k, it will not have any effect on the Wigner quasi-distribution function neither will it have an 
effect on the measurable quantities. This is follows from the cancellation of the additional phase factors in equation 
(f79|). Hence in order to cause a change in the physical system, some other factors should be changed. This statement 
applies also to the Bloch representation. To clarify, we give the connection between the fcg-representation and the 
Bloch representation, 

a{k,q)=Y J B n {k)^ nk {q) (84) 

n 

where the coefficients B n (k) are the wave functions in the Bloch representation, and ^ n k(q) are the Bloch functions. If 
a change of the quasi-momentum is necessary which will also change the momentum one should introduce a function of 
the quasi-coordinate and the quasi-momentum, f(k, q) in the definition of a(k, q). Usually the function exp J A(fc, q)dt 
is used with A(fc, q) chosen to be the dispersion relation of the electrons with the homogeneous fields set to zero. The 
issue of the validity of this procedure is reviewed by Nenciu (1991), and we shall not dwell on it here. All that has 
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been said so far, is that a change in the quasi-momentum cannot by itself lead to a kinetic change in the physical 
system. Hence we suggest that a change in the quasi-momentum due to external field is a gauge effect. 

It is sometimes stated, that the dynamics of an electron in a lattice can be described by the so called acceleration 
theorem which is generally written in the form 

dk 

— = -eE - {e/c)v av x H (85) 
at 

where k is the quasi-momentum of the particle and v av is the average velocity of the electron. It is usually stated 
that one may find a basis of eigenfunctions in such a way that it will make the average velocity of the electron, after 
a sufficient long time, to be v a „ = V fc A(k), where A(k) is the dispersion relation (Adams and Argyres, 1956; Wilson, 
1965). 

The acceleration theorem is regarded as a quasi-classical approximation that can be used only if the external fields 
are weak and if the movement of the Bloch electron is restricted to a small part of the Brillouin zone. The existing 
derivations of the acceleration theorem (Wilson, 1965; Landau and Lifshitz, 1984; Kittel, 1963; Kohn, 1959) are very 
complicated, controversial and are not rigorous. Hence, their validity is not well defined. An exception is the paper 
of Zak (1972) who has derived the acceleration theorem for an electric field in a rigorous way. 

It is interesting to note a similarity between the acceleration theorem and equation (p7j). Here v Q1I corresponds to 
the velocity of the observer da;"/dt" and the canonical momentum to the quasi-momentum. To make it more clear 

we remark that P can (X; x") is equal to the expectation value of the momentum p = k + K. Remembering that the 
expectation value of K can take only discrete values, we realize that there are short periods of time At" , when the 
change of the expectation value of the momentum is equal to the change of the expectation value of k. At these 
periods of time the acceleration theorem and equation (|27j) are identical. 



C. Canonical energy change of a heavy mass polaron in electric field 

The problem of a heavy mass polaron^] (van Haeringen W, 1965) is solved with the help of the quantum mechanical 
equivalent of the classical Hamilton- Jacobi equation. We shall solve here the problem of a heavy mass polaron in a 
constant and homogeneous electric field. This example is also interesting in and of itself, since this problem has never 
been solved by means of a Dyson equation, (see discussion in the book of Mahan (1990), p. 509-510). 

We are concerned with a homogeneous system with a single electron in a conduction band and assume an electron- 
phonon interaction with optical phonons. For that kind of system G- — and only one Green's function should 
be evaluated (in the general case there are two independent Green's functions) and we shall take it to be G c . This 
function satisfies the quantum Hamilton- Jacobi equation (pOl) that for our system takes the form 



1 

2~7rft 



£-e p + ^-(V x + eEd s ) 2 
8m 



G c (p, £ ,X) 



l + E c (p,e,X)G c (p,e,X) (86) 



with £ p = p 2 /(2m). The right hand side of the equation takes this simple form since the polaron has a heavy mass 
and all the derivatives of the Green's function with respect to the momenta are assumed to vanish. Hence, the phase 
loop function (that represents kinetic effects) vanishes. The self-energy term has the general form 

S c (p, e)= J dc^dqG c (p, e-u- £ q )£ c (q, w)r(p, e, q, w) (87) 

D c (q, u>) is the phonon Green's function and the vertex part is given by the Ward identity (Engelsberg and Schrieffer, 
1963) 

r( P ,e,q,o,) - -— 



3 The heavy mass polaron is discussed only for the sake of a demonstration. We therefore do not provide an introduction to 
the problem of polarons and assume that the reader is familiar with the problem. The background can be found in text-books 
(see, e.g., Mahan, 1990). 
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Fourier transforming equation (pq), one obtains 

d 



-ih- - s p 



— (-ihV x + eEt) - 

8m v / 



G c (p,i) = hS(t). 



(89) 



where 



$(i)= J dTdq£> c (q, T)e fe i T 

The solution of this equation is 

G c (p,t,X) =0(i)exp i 



(e p + eEX)i+ / dr*(r)) 



(90) 



As can be seen from this equation the energy of the polaron is composed of three parts corresponding to its kinetic, 
potential and phonon interaction energies. It is instructive to search for a solution of equation (|8^) that does not 
depend on X 



1 



8m 



G c (p,e) = l + S c (p,e)G c (p )£ ). 



It reads 



G c (p,i) =0(t)exp 



T £ot - 



24m 



t° + 



dT$(T) 



(91) 



(92) 



Going back to equation (27) we see that the change of the canonical momentum of the particle at rest in a static and 
homogeneous electric field is eEt. The change of the canonical energy due to the electric field is 



&can (^) 



/ dt ^ Ep ( f ')=s(I E )7' <lt ' ti 4(I E ) , f 



Integrating over time t we arrive at the phase factor in equation (Ml) which can be now written in the form 



G c (p,i) = 0(t)exp 



UJ (e p -e can (T) + ^(T))dT 



(93) 



It does not contain the potential energy term. There is a new term instead which corresponds to the energy that a 
free particle would acquire if it were moving freely in the electric field. 



VI. SUMMARY 



In this paper we formulate a gauge-invariant formalism for transport processes starting from the definition of the 
gauge-invariant Wigner quasi-distribution function, we also derive quantum transport equations and their hydro- 
dynamic representation. Throughout the paper we emphasized the kinetic vs canonical meaning of the quantities 
discussed. We believe that this formulation will both simplify studying transport phenomena and reveal new as- 
pects that are not easily shown with other available formalisms (e.g., quasi-classical Boltzmann equation or Kubo 
correlators). 



APPENDIX A 



The equality 




2G 



1 



ic: 



exp | — / y'A^X + sy)ds | (A.l) 



in which the canonical and kinetic momenta are connected by the equation 



Pkin = Pcan - ~A{X) = -iTl— - -A(X). (A.2) 

c oX c 



is derived. 

First, the operator 



U(s) = exp L'-sP^y^ exp UsP&J^ (A.3) 

is defined. It satisfies the following two conditions 

U(s)\ s=0 = 1 (A.4) 



and 

dU(s) ie 
ds ch 

where 



rM* + sy)U(s) (A.5) 



A(X + sy) = exp ( --sP^J^ j A(X) exp f ) . (A.6) 



Then integrating equation ( A.E) the equality (A.l) is obtained 



APPENDIX B 



The phase loop function was introduced in our previous paper (Levanda and Fleurov, 1994). Here we consider its 
shape for systems in a static and homogeneous field. The phase loop function is another representation of the slashed 
derivatives appearing in equations ( |30| ) and (|3l]) and discussed in detail in main body of the paper. 

The phase loop function corresponding to an electron loop with three vertices is defined by the equation 



M(K\3,F(X))=M({m}\F(X)) 



J d 4 Xl d 4 x 2 d 4 x 3 exp f I tfxiA ■ exp (jgb ({a*})) (B.l) 



i=l 

where 

3 X ' 2 



$({£fc}) = X>fc / A^(X k + Sk y k )ds k , (B.2) 

i=1 -1/2 

y k = x k+ \ — x k , X k — (xfe+i + x k )/2, and cyclic conditions x^+i = Xi being assumed. Such a phase loop function 
depends on three four-dimensional variables {ir^ }f =1 , on the electromagnetic field tensor F XfJl (X) and is explicitly 
gauge invariant. In the absence of electromagnetic fields the phase loop function takes its trivial form 

M^\3,F(X)) = ]JS(4), (B.3) 
Choosing x" — xi in equation ( |l7| ) we find that 
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$ ({Xk}) = [{x 2 ~ Xl)n(X3 ~ Xl)\ ~ {x 3 ~ Xl)n(X2 - Xl)\] . (B.4) 

This is just an expression for the flux through a triangle that has its three vertices at the points 11,0:2,2:3, 

$({x fc })= J F^dXpAdxx, (B.5) 

S{xi,X 2 ,X 3 } 

where S{x\, x%, £3} is the surface covering this triangle. 

The integration in the definition of the phase loop function is now straightforward and results in 

M(tt I 3,^ A ) = 

S (j2 WM ex P [(-ieh/4c)F li xff r " t a*> x ] n Sfa/h). (B.6) 

\i=l J i=2 



The latter equation (B.6) can be easily generalized to a N vertex electron loop in a homogeneous static field. The 
polygon with N vertices can be subdivided into N — 2 triangles and the phase loop function for such N vertex electron 
loop is 

M(tt I N,^ x ) = 




exp 



AT-l 



k=2 



N 



HSinJh). (B.7) 



i=2 



Comparing this equation with equation (B.3), one realizes that the momentum or energy are no longer conserved 
at each vertex. In this representation, the phase loop function shows in a simple way, how the field interferes with the 
four-momentum conservation. Derivatives of the delta function mean that derivatives of the Green's functions and of 
the vertex functions appear in analytical expressions for the diagrams. These derivatives stand for the electron recoil 



due to the field (compare this to the recoil expansion of van Hearingen, 1965). Equation (B.7) converges to equation 



(B.3) as J~^\ approaches zero. 
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